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Abstract 

We have studied how parton distributions based on the inclusion of nonlinear scale 
evolution and constraints from HERA data affect charm production in pp collisions 
at center-of-mass energies of 5.5, 8.8 and 14 TeV. We find that, while the resulting 
enhancement can be substantial, it is very sensitive to the charm quark mass and the 
scale entering the parton densities and the strong coupling constant. 
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1 Introduction 



Global fits of the parton distribution functions (PDFs) such as those by CTEQ [H E] 
and MRST 0HIH1, based on the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 

scale evolution, successfully describe the proton structure function, F 2 (x,Q 2 ), deep 
inelastic scattering (DIS) data in the "high (Q 2 , x)" region, Q 2 > 10 GeV 2 and x > 0.005. 
However, it is not possible to maintain the excellent high (Q 2 ,x) DIS fit while si- 
multaneously fitting the "low (Q 2 ,x)" region, 1.5<Q 2 <10 GeV 2 and x< 0.005 In 
addition, the next-to-leading order (NLO) gluon distribution becomes negative for suf- 
ficiently small x at the few GeV 2 scales. 

Nonlinear corrections to the PDF evolution based on gluon recombination were first 
derived by Gribov, Levin and Ryskin [Zj as well as Mueller and Qiu [H]. Recent work 
PI showed that adding these GLRMQ terms to the DGLAP equations can improve the 
overall leading order (LO) fits to the HERA DIS data [1 0J . The rapid Q 2 evolution in 
the low (Q 2 ,x) region from DGLAP alone is slowed by the GLRMQ recombinations. 
At Q 2 scales far above the initial scale Ql, the Q 2 evolution of the PDFs is again 
described by the DGLAP equations since the GLRMQ terms become negligible, see 

Fig.m 

While the quark distributions are directly constrained by the HERA F 2 (x, Q 2 ) data, 
the gluon distribution is constrained by the F 2 slope, dF 2 /d\iaQ 2 . The fairly modest 
measured slopes force the LO DGLAP gluon distributions to be nearly independent of 
x for scales of a few GeV 2 [HOIl and also force the NLO gluons to be negative. When the 
nonlinear terms are included, the slowing of the Q 2 evolution leads to an enhancement 
of the small-x gluon distributions at Q 2 ^ 10 GeV 2 relative to the LO DGLAP gluon 
distributions, subject to the same constraints from HERA Pj. The effect of nonlinear 
evolution on the NLO distributions is not yet fully explored by global fits which include 
the low (Q 2 , x) region but the results of Ref. (Hj suggest that the enhancement is smaller 
than at LO. This effect alone thus seems unlikely to produce positive NLO small x gluon 
distributions at the few GeV 2 scales. 

In spite of the problems described above, the quality of the global DGLAP fits to 
the HERA data [H 0] is good. The y 2 per degree of freedom is close to one even when 
the low (Q 2 ,x) region is included. Therefore, F 2 measurements at HERA alone may 
not clearly differentiate DGLAP from nonlinear evolution and more direct probes of the 
gluon distribution are needed. In this paper, we study whether the parton distribution 
functions generated with the LO DGLAP+GLRMQ evolution in Ref. |Hj could give rise 
to any significant enhancement relative to the DGLAP-evolved PDFs in charm quark 
hadroproduction. Charm production is the best candidate process since the charm 
quark mass is relatively low, 1.2<m c <1.8 GeV, and its production is dominated by 
gluons. These two characteristics should lead to the most favorable conditions for 
a possible effect. The Q 2 scale at which the total cross sections are calculated is 
proportional to m 2 and 4m 2 . Thus the results are sensitive to Q 2 . Unfortunately, due 
to the small charm mass, the scale dependence is still significant at NLO [TT| . 

We focus on pp collisions since these nonlinear distributions are not yet available for 
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nuclei. Our calculations are at leading order only since the nonlinear parton distribution 
functions, referred to as EHKQS hereafter, are only evolved to LO in Ref. [0]. We 
calculate the possible effect as a function of rapidity and transverse momentum of 
the charm quark and the cc pair invariant mass. We study collisions at center-of-mass 
energies VS = 5.5, 8.8 and 14 TeV at the Large Hadron Collider (LHC). These energies 
correspond to the planned per nucleon energies of Pb+Pb, pPb and pp collisions, 
effectively spanning the LHC energy regime. 

2 Formalism and inputs 

2.1 Cross sections and parton distribution functions 

Inclusive differential charm cross sections at high energies are, to first approximation, 
computable assuming factorization. The cross section may be expressed as 

d(Tpp^ C cx(Q 2 , y^) = fi( x i,Q 2 ) ® fj(x 2 ,Q 2 ) ® da ij ^c 5k (Q 2 ,x 1 ,x 2 ) (1) 

i,j,k=q,q,g 

where o-ij^ C ck(Q 2 ,Xi,x 2 ) are the perturbatively calculable partonic cross sections for 
charm production at scales Q 2 ^> Aq CD , x 1 and x 2 are the momentum fractions of 
the partons involved in the hard scattering and fi(x,Q 2 ) are the free proton parton 
densities. At LO, only the gg and qq channels are available and k = 0, i.e. no other 
partons are produced with the cc. The LO matrix elements and partonic cross sections 
for charm production can be found in Ref. ^2]. In the following, we consider the 
triple differential distributions, d 3 a/(dp T dydy 2 ), where y and y 2 are the rapidities of 
the quark and antiquark and pr is the quark transverse momentum. We also study the 
inclusive distributions da/dy, da/dpr and da/dM where M 2 = 2m^(l + cosh(y — y 2 )) 
is the square of the pair invariant mass and m\ = p\ + m 2 c is the transverse mass of 
the quark. 

The new inputs in this straightforward calculation are the nonlinearly-evolved pro- 
ton PDFs. The EHKQS sets 4 in Ref. |9] employ the CTEQ5L [TJ and CTEQ6L P 
PDFs as a baseline and require a good fit to the HERA F$(x, Q 2 ) data over the full 
(Q 2 ,x) range Three EHKQS sets were obtained when the GLRMQ terms were 
included in the analysis. All three EHKQS sets employ the initial scale Qq = 1.4 
GeV 2 and Aq CD = 0.192 GeV for four flavors but differ in the treatment of the charm 
mass threshold in the evolution. In set 1, a nonzero charm distribution was allowed 
at Qq. Sets 2a and 2b assumed that there was no charm quark distribution below Qq, 
turning on when Q 2 = m 2 c > Qq. Set 2a assumed m c = 1.3 GeV while set 2b took 
m c = Qq = a/L4 GeV. The input gluon distribution is the same in all cases so that 
the small differences at high Q 2 arise from the treatment of the charm quark evolution. 
The choice of set 1, set 2a or set 2b therefore makes very little difference in the overall 
effect of the nonlinear terms on charm production. Thus we only use EHKQS set 1. 

4 available at www.urhic.phys.jyu.fi 
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We work at leading order since the EHKQS sets are only evolved to LO using a 
one-loop evaluation of the strong coupling constant a s . Thus these LO distributions 
should generally not be mixed with NLO matrix elements and the two-loop a s . The 
charm quark pt distribution is broadened at NLO relative to the LO calculation [TT] . 
Therefore we study the ratios of calculations with EHKQS relative to those with a 
standard LO PDF set evolved using the DGLAP equations alone. We quantify the 
effect with respect to the CTEQ61L parameterization, the most recent LO fit to the 
PDFs that also uses a one-loop a s PJ. The minimum scale of CTEQ61L is Ql = 1.69 
GeV 2 . This LO fit obtained a slightly higher value of Aqcd for four flavors, 0.215 MeV. 
In our CTEQ61L calculations, for consistency with the PDF set, we use this value in 



2.2 Comparison of EHKQS and CTEQ61L gluon distributions 

Before presenting our results, it is instructive to discuss the differences between the 
EHKQS and CTEQ61L gluon distributions in more detail. Since the high \/S collisions 
studied here probe the very low x region, some remarks on the region of applicability 
of the PDFs are in order. Below the minimum x and Q 2 values assumed in the fits, the 
PDFs are essentially unconstrained. At the scales studied here, described in the follow- 
ing section, Q 2 always remains above the minimum scale of the PDF sets. However, 
the region below the minimum x, x mm , is reached at large rapidities and intermediate 
scales at the LHC. Thus the behavior of the PDFs below x min is an uncertainty for all 
PDFs. This situation can only be improved by more extensive small x constraints on 
the PDFs. 

The minimum x of the EHKQS sets is a;™ KQS = 10" 5 . For x < s™ KQS and Q 2 of 
a few GeV 2 , neglected power-suppressed terms in the evolution become important and 
the DGLAP+GLRMQ results are no longer trustworthy p]. The CTEQ61L minimum 
is an order of magnitude smaller than £™ KQS , ^ m in EQ61L = 10~ 6 . The very small x 
regions below x m \ n are not excluded from our calculations. Instead, we assume that 
below x mm , fi(x < x mm , Q 2 ) = fi(x m i n , Q 2 ) for each set. We note that in the CTEQ61L 
table [TJ] the distributions are not constant below x m ^ EQ61L . 

It is illustrative to compare the EHKQS set 1 and CTEQ61L gluon distributions 
as a function of Q 2 for several values of x, shown in Fig. [TJ Due to the nonlinear 
evolution, the Q 2 dependence of EHKQS set 1 is rather mild compared to CTEQ61L 
which approaches a constant at Qq and x — > 0. Note that for all x > 1CT 5 , the EHKQS 
distributions are always above CTEQ61L although, at higher scales, the distributions 
lie very close together. In the unconstrained region where x < 10~ 5 , the situation 
clearly depends on how the extrapolation towards x — > is done. The CTEQ61L 
parameterization continues to x m ™^ 61L . In this very low x region, at Q 2 = 4 GeV 2 , 
the CTEQ61L gluon distribution at the EHKQS distribution, fixed 

at the value of £™ KQS , and continues to rise. Therefore, the behavior of the relative 
kinematics distributions we compute can be very sensitive to the treatment of the 
unconstrained x region. Since the two distributions are also quite sensitive to the 
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Figure 1: Comparison of the EHKQS set 1 (solid curves) and CTEQ61L (dashed curves) 
gluon distributions as a function of Q 2 for, from lowest to highest, x = 10 -2 , 10 -3 , 10 , 
1(T 5 and, for CTEQ61L only, 1(T 6 . 

scale, the ratios we compute will also be strongly scale dependent. 
2.3 Scale choice 

In our calculations, we use values of the charm quark mass and scale that have been 
fit to the total cross section data using NLO calculations. The total cross section data 
cannot be fit by adjusting m c and Q 2 with a full LO calculation, employing LO PDFs 
and the one-loop ot s , because the resulting m c would be too small for perturbative 
applications. See Ref. (TTj for more discussion. The best agreement with the total 
cross section data is obtained with m c = 1.2 GeV and Q 2 = Am 2 for standard DGLAP- 
evolved NLO PDFs such as CTEQ6M [2| and MRST Nearly equivalent agreement 
may be obtained with m c = 1.3 GeV and Q 2 = m 2 (TU [TJj. Thus our main results 
are based on these inputs. Alternative fits to the fixed-target total cross sections can 
be achieved with larger values of m c by separating the factorization scale, Qf, from 
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the renormalization scale, Qr, e.g. Q 2 F = 4m 2 and Q 2 R = m 2 c or m 2 / 4, and allowing 
the fast running of a s in this Q 2 R region to increase the cross sections. Note that if 
Q 2 = Q 2 p < Qq, the PDFs are unconstrained in Q 2 . We keep Q 2 F = Q 2 R , as in all 
typical PDF fits such as Refs. [21 CHI, limiting ourselves to relatively small values of m c 
to obtain agreement with the total cross section data. The PDFs are thus evaluated 
above Qq for the masses and scales we use. 

Note that we have discussed scales proportional to m 2 in the calculations of the 
total cross sections. Such scales are used because the total NLO partonic cross section 
can be written analytically as a function of Am 2 /s where s is the partonic center-of- 
mass energy squared |18| . In this case, the charm quark mass in the only relevant 
scale. However, in inclusive distributions such as we compute here, the quark px also 
enters since a scale proportional to m\ is needed to control pr-dependent logarithms. 
Therefore, in our calculated ratios of distributions, we take Q 2 = with m c = 1.3 
GeV and 4m| with m c = 1.2 GeV. 

Whether the high and low energy behavior of the charm cross section can be de- 
scribed simultaneously by the same values of m c and Q 2 is an open question. Since 
the slopes of the gluon distributions differ at low and high x, the scale dependence is a 
strong function of v^S*- At high x, xf g (x, Q 2 ) is larger at low Q 2 than at higher scales. 
Thus at fixed target energies, aim 2 ) > a(4m 2 ). At collider energies, such as at the 
LHC, x is small and xf g (x,Q 2 ) is increasing with Q 2 so that <r(4m^) > cr{m 2 ) even 
though a s {m 2 ) > a s (4m 2 ). We thus extend the parameter space of our calculations to 
study the possible effect on higher masses, m c = 1.8 GeV with Q 2 = m F and 4m F . 



3 Results 

We start with results insensitive to the unconstrained region at x < x min = 1CT 5 . 
In inclusive kinematics with an identified charm quark and fixed xt = 1m F j^S^ the 
unconstrained x-region contributes to charm production at high rapidities, in the region 

y u = In (l/x T - \jl/x 2 T - l/zmin^ < \y\ < In (l/x T + \JlJx\ - l/aw) • (2) 

The upper limit, close to the phase space boundary, is not of interest here. Expanding 
the lower limit, y u , in powers of x^/x min 1, we arrive at y u ~ \n[m T / (x min ^/S)] > 
ln[m c /(x min V / 5')]- If rn c = 1.2 GeV, the small x region contributes to charm produc- 
tion at \y\~y u = 2.2, 2.6 and 3.1 for = 14, 8.8 and 5.5 TeV, respectively. If 
m c = 1.8 GeV, the unconstrained region is probed when \y\~y u = 2.6, 3.0 and 3.5, 
respectively. 

First, we study the ratio of the fully differential cross sections calculated with the 
EHKQS densities relative to the CTEQ61L densities at fixed y and y2, 

R{pt, y, V2) = vrn a a \ 3 

d 6 a {CTEQ61LJ / \dpTdydy2) 
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Figure 2: We present R(pT,y,U2) for fixed y and j/2 as a function of charm quark px at 
\/S = 14 TeV, (a) and (d), 8.8 TeV, (b) and (e), and 5.5 TeV, (c) and (f), in pp collisions. 
The results are shown for y = y2 = (solid curves), y = 2 y 2 = (dashed) and y = y 2 = 2 
(dot-dashed). We show m c = 1.2 GeV and Q 2 = 4m 2 - on the left-hand side and m c = 1.3 GeV 
and Q 2 = m 2 - on the right-hand side. Note the different scales on the left- and right-hand 
axes. 



In Figs. [2] and El we present this unintegrated ratio as a function of pr for several values 
of y and y 2 , all within the range constrained by data: y = y 2 = (solid curves), y — 2, 
2/2 = (dashed) and y — y 2 — 2 (dot-dashed). Figure El presents i?(pr, 2/, 2/2) for the 
values of mass and scale that agree best with the total cross section data, m c = 1.2 
GeV and Q 2 = 4ni^ on the left-hand side and m c = 1.3 GeV, Q 2 = m 2 ^ on the right- 
hand side. Figure El shows the same unintegrated ratios for our upper limit on m c , 1.8 
GeV and Q 2 = m^^Am 2 ^. The center-of-mass energies shown are v^S" = 14 (upper), 8.8 
(middle) and 5.5 (lower) TeV. 

Figures [2] and El which demonstrate the sensitivity of the enhancement to m c and 
Q 2 , can be understood by inspection of Fig. [TJ At these high energies, the gg channel 
dominates so that the partonic cross sections essentially drop out of the ratios, and 

r> ( \^ x lfg( x l>Q 2 )ERKQS X 2 f g (x 2l Q 2 )eHKQS _ d / n .2\ r> I rf2\ , A \ 

R(pt, y, V2) ~ - 7- 2) , 2] = Rg(Xl, Q )R g (x 2 , Q ) , (4) 

^l/g^ljV JCTEQ61L x 2J g { x 2: V JCTEQ61L 

where we have denoted the ratio of the EHQKS and CTEQ61L gluons by R g . The x val- 
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Figure 3: As in Fig.EJbut for m c = 1.8 GeV and Q = Am\ (left-hand side) and for m c = 1.8 
GeV and Q 2 = m\ (right-hand side). Note the different scales on the left- and right-hand 
axes. 



ues are easily calculated at LO for each curve using the definitions x\p = rriT[exp(±y) + 
exp(±y 2 )} / VS . Since decreasing a/5" increases x and consequently decreases R g {x, Q 2 ), 
the enhancement decreases with energy. In addition, R g (x,Q 2 ) decreases with increas- 
ing Q 2 , so that increasing m c and Q 2 also reduces the enhancement. Both x\ and x 2 
are small when y = y 2 = so that the enhancement is largest at midrapidity. Moving 
away from midrapidity, e.g. to y, y 2 > 0, increases x\ and decreases x 2 , correspondingly 
decreasing R g (xi,Q 2 ) and increasing R g (x 2 ,Q 2 ). The rapidity dependence shown on 
the right-hand sides of Figs. [2] and El for the lower Q 2 , is rather weak. The difference 
between the dashed (y = 2,y 2 = 0) and dot-dashed (y = y 2 = 2) curves is very small 
and the enhancement at y = y 2 = 2 lies marginally above that for y = 2, y 2 = over 
all pt- When the scale is small, the CTEQ61L gluon distribution changes very slowly 
with x for x < 0.01, about 20% for 10~ 4 < x < 10~ 3 when Q 2 = ml and m c = 1.3 
GeV (pt ~ 0), so that R(pT,y,y 2 ) is not a strong function of y and y 2 . However, 
the results for the larger scales, shown on the left-hand sides of Figs. [2] and exhibit 
the opposite behavior along with a stronger rapidity dependence. At larger scales, the 
slope of the CTEQ61L gluon distribution with x is considerably stronger, resulting in 
a factor of two difference in the gluon distribution over the range 10~ 4 < x < 10~ 3 
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Figure 4: We present R(y), (a) and (d), R(pr), (b) and (e), and R(M), (c) and (f), in pp 
collisions at y/S = 14 (solid), 8.8 (dashed) and 5.5 (dot-dashed) TeV. The left-hand side shows 
m c = 1.2 GeV and Q 2 = 4m^, the right-hand side m c = 1.3 GeV and Q 2 = m T . 



when Q 2 = 4m 2 and m c = 1.2 GeV (pr ~ 0), introducing a stronger rapidity depen- 
dence of R[[pr, V, 2/2) at higher scales. The nonlinearities die out at large scales since 
the EHKQS gluons become similar to the CTEQ61L gluons so that the enhancement 
in R(pT,y,y2) disappears at large px- The ratio does not become equal to 1 because 
CTEQ61L and the EHKQS sets are not identical either at high Q 2 or larger x. Note 
that a^(EHKQs)/«g(cTEQ6iL) w 0.9, allowing R(pT,y,V2) to drop below 1 at high px- 
Next, we turn to the integrated ratios, 

_ dcr (ebkqs) / dy ^, . _ do-(EnKQs)/dp T R(M) ~ < ^ <T ( EHKQS )/^ J ^ 

da(cTEQ61L)/(i?/' ' ' ' C?Cr(cTEQ61L)/(ij9T' da (cTEQGlh) / dM ' 

(5) 

shown in Figs. 0] and El All three energies, y/S = 14 (solid), 8.8 (dashed) and 5.5 
(dot-dashed) TeV, are shown in each plot. Figure 0] shows the integrated ratios for the 
mass and scale values that best agree with the total cross section data, m c = 1.2 GeV, 
Q 2 = Am? on the left-hand side and m c = 1.3 GeV, Q 2 = m T on the right-hand side. 
Figure El shows the corresponding results for m c = 1.8 GeV with Q 2 = Am? (left-hand 
side) and (right-hand side). 



8 




5.0 -5.0 -2.5 



0.0 

y 



2.5 5.0 



1 

\y 


1 1 1 1 

m c = 1.8 


l i 

GeV - 




Q 2 = m T 2 






. (e) 





— 1.6 



— 1.4 



1.2 
1.0 



10.0 0.0 2.5 



5.0 

P T (GeV) 



7.5 10.0 



5 M (GeV) 10 




5 M (GeV) 10 



Figure 5: As in Fig. Q] but for m c = 1.8 GeV and Q 2 = 4m 2 - (left-hand side) and Q 2 
(right-hand side). 
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As discussed previously, at midrapidity the results for R(y) are insensitive to the 
EHKQS extrapolation region x < x™ KQS . The magnitude of the enhancement in R(y) 
can be understood from the unintegrated results. Since R(y) is integrated over p?, 
it does not only reflect the enhancement at % = m c (p T = 0) because Q 2 oc m\ 
[TO] and the pt distribution peaks around px ~ 1 GeV. When R(y) is calculated with 
Q 2 = m?, shown on the right-hand sides of the figures, the ratios are broad due to the 
stronger growth of R g (x2, Q 2 ) relative to the reduction of R g (xx, Q 2 ) with increasing y, 
i.e. decreasing xi and increasing x\. The ratio is broad because the CTEQ61L gluon 
distribution is relatively flat as a function of x for Q 2 ~ 2 — 3 GeV 2 . The enhancement 
decreases and broadens with decreasing energy. 

The rather sharp turnover in R(y) indicates where the extrapolation region x < 1CT 5 
begins to contribute. The rapidity at which x < 10~ 5 is larger for lower energies and 
larger m c , see Eq. (J2J). The clear decrease of R(y) below 1 at large rapidity for Q 2 = Am% 
is perhaps surprising. This effect is a consequence of the small-x extrapolation adopted 
for the PDFs and can be understood through an examination of Fig. [TJ While the 
EHKQS gluon distribution is fixed at its value at x™ KQS = 1(T 5 , the CTEQ61L 
distribution continues to change until x|^™ Q61L = 10~ 6 . At v^S" = 14 TeV, the rapidity 
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at which some x values fall below x min is y u = 2.2 (2.6) for m c = 1.2 (1.8) GeV. For 
\y\ > y u and Q 2 > 4 GeV 2 , as x decreases, the CTEQ61L gluon distribution increases 
considerably above the EHKQS distribution. Thus R{y) < 1 at large rapidities. 

Since the rapidity distributions are rather flat, there are still important contribu- 
tions to the pt and mass distributions, up to ~ 30% and 40% respectively at y/S = 14 
TeV for m c = 1.2 and Q 2 = "im 2 , from the extrapolation region. The contribution 
from the extrapolation region to R(M) is larger since when m\ is small, M can still 
be large because the difference y — y<i can be large while either or both y and y<i can 
be in the low x region. Thus the sensitivity of R{pt) and R(M) to the unconstrained 
region should be kept in mind. 

The magnitudes of R{j>t) and R(M) in Figs. |U and El are similar to those of 
R(Pt, V, Vi) in Figs. [2] and El However, the relative energy dependence is reversed 
from the large to small scales. At the largest v^, the contribution from the x < 10~ 5 
region is the largest, and as seen in Fig. [TJ if the scales are e.g. above 5 GeV 2 , the 
CTEQ61L gluon density becomes higher than the EHKQS gluon density at x < 10~ 5 , 
suppressing R(pr) relative to R(pT,y,y2)- At smaller y/S the higher x values reduce 
the contribution to da/dpT from the region x < 10~ 5 . Thus the suppression of R{pr) 
from the extrapolation region is reduced with decreasing energy, explaining the relative 
values of R(pt) and R(M) with energy shown with Q 2 = 4m^ (the left-hand sides of 
Figs. |U and EJ). In contrast, for Q 2 = m 2 -, the EHKQS gluon density is always higher 
than CTEQ61L at small x (see Fig. [l]), and the energy dependence remains the same 
as in Figs. [2] and 0J 



4 Discussion 

Including the GLRMQ terms with DGLAP evolution enhances the low x proton gluon 
distribution jHj. It may, however, be very difficult to distinguish between linear and 
nonlinear Q 2 evolution on the basis of the F 2 HERA data alone. Other probes of 
the low x gluon distribution, such as charm production in DIS at HERA and in pp 
collisions at the LHC, will hopefully provide the necessary further constraints. In 
this paper, we have demonstrated how the nonlinear PDFs constrained by the HERA 
data can cause a significant enhancement in the LO charm quark cross sections in pp 
collisions at the LHC. The enhancement is defined relative to the results expected with 
the pure DGLAP PDFs which also fit the same HERA data. Quantitatively, however, 
this enhancement is shown to be very sensitive to the charm quark mass and the scale. 
Clearly, in addition to collecting more data, further theoretical input, such as high 
energy resummation of the heavy quark cross sections j20j, is needed to reduce this 
sensitivity. 

Our basic message is the following. Currently the high-precision HERA data on 
the structure function F 2 (x, Q 2 ) and its derivative with respect to Q 2 provide quite 
stringent constraints on the gluon distribution. Consequently, the CTEQ61L and 
MRST2001LO gluon distributions are very similar. If a significant enhancement of 
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charm production, unexpected from DGLAP evolution alone, is found at the LHC, 
such an enhancement cannot be absorbed into the DGLAP-evolved gluon distributions 
without introducing a steeper InQ 2 slope of F 2 (x,Q 2 ) than allowed by the HERA 
data. Therefore, such an enhancement would be a signal of nonlinear effects on the 
PDF evolution. 

The unintegrated ratios, R(pT,y,y2), at \y,U2\ < 2 and pr ~ are quite large for 
the masses and scales that best agree with the fixed target data. As discussed above, 
in charm production at central rapidities the PDFs are constrained by the HERA data. 
For the smaller charm mass and larger scale, m c = 1.2 GeV and Q 2 = 4m^, we find 
R(pT,y,y2) ~ 1.4— 1.5. The enhancement is even larger for the smaller Q 2 , ~ 3.7 — 5.5 
with m c = 1.3 GeV and Q 2 = m T . The corresponding enhancement factors are smaller 
for the upper limit on the charm mass, m c = 1.8 GeV, 1.1 — 1.2 for Q 2 = 4m 2 - and 
1.7 — 2.1 for Q 2 = m 2 -. We have also shown that the enhancement disappears with 
increasing pt due to the decreasing importance of the nonlinear terms at larger scales 
in the PDF Q 2 evolution. 

The enhancements of the integrated ratios are somewhat reduced from the uninte- 
grated results above. The rapidity enhancement is R(y) ~ 1.07 — 1.8 at midrapidity 
for the energies studied. For the pt enhancement, we find R(pr = 0) ~ 1.05 — 1.25 
with Q 2 = 4m T and ~ 1.6 — 4.5 with Q 2 = over the range of m c we investigate. 
The mass dependent enhancement, R(M = 2m c ), is similar albeit a bit smaller than 
R(Pt = 0). At larger px and M, the effect dies out rather quickly, see Figs. 0] and 03 
The main uncertainty in R(pr) arid R(M) is the sizable contribution from the region 
x < 10~ 5 , currently unconstrained by the HERA data. Thus the enhancement also 
depends on the extrapolation of the PDFs in this x region. 

The enhancement we calculate here should be an upper limit on the possible ef- 
fect at LO, as we now discuss. The recombination radius of the gluon ladders in the 
GLRMQ terms was assumed to be the proton radius in the EHKQS analysis The 
proton radius was the lower limit on the recombination radius since it gave the strongest 
possible recombination effect while still describing the HI data [10] over the full range, 
x > 3 x 10~ 5 and Q 2 > 1.5 GeV 2 , without entering the saturation region. In the sat- 
uration region, further nonlinear terms become important and the DGLAP+GLRMQ 
evolution breaks down. Thus, in this approximation, the EHKQS gluon distribution is 
the upper limit on the LO gluon distribution. The results obtained here are then an 
upper limit on the enhancement. Studies of the nonlinear PDFs, particularly in the 
context of NLO DGLAP evolution, are needed to go beyond this approximation. 

Finally, we discuss possible detection of this enhancement. Fragmentation and 
decay should not wipe out the effect, as also seen for shadowing in pA collisions [21]. 
Although the enhancement R(pT,y,y2) is largest, see Figs. [2] and 01 dileptons may not 
be the best channel to measure the enhancement because the origin of the individual 
pairs is unknown. The lepton px and y also does not correspond directly to the quark 
Pt and y. In addition, it is not clear how the enhancement survives a minimum pt cut 
of the charm decay leptons, even though these decay leptons are a major component 
of the dilepton spectrum. Further simulations of this would be worthwhile. Since 
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the enhancement disappears at large px, the smallest possible pr cut is desirable. 
At higher lepton pair masses, bb decays dominate the dilepton mass distributions. 
The enhancement rapidly decreases with increasing quark mass so that effects on bb 
production are very small. 

Single leptons from charm decays [22] may be a better possibility. The single lepton 
enhancement may, however, be somewhat reduced relative to that of the pair. Of course 
fully reconstructed D mesons would be the most desirable option. Reconstruction 
should be possible for pt ~ 1 GeV or less in ALICE with identified kaons - More 
complete simulations with nonlinear PDFs should thus be performed. 

Ideally, a better place to search for the nonlinear effects would be pA collisions at 
the LHC |21j. The GLRMQ corrections are expected to get an oc A 1 ^ enhancement 
from the nuclear size, causing the nonlinearities to be significant at larger values of x 
and Q 2 than in free protons, see Ref. [2S1 in Ref. |24j. Thus the enhancement should 
be more pronounced for nuclei. Before computing charm cross sections in nuclear 
collisions, however, the nuclear PDFs should be analyzed within the DGLAP+GLRMQ 
framework, including constraints from nuclear deep inelastic scattering. 
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